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We study the effect of contact hyperfine interaction on the nuclear spin diffusion coefficients in 
semiconductor quantum dots. The diffusion coefficients are calculated with both the method of 
moment and density matrix. We show that nuclear spin diffusion is strongly suppressed by the 
nonuniform hyperfine coupling resulting from the confined electron wavefunction. Our results agree 
with the observed suppression of nuclear spin diffusion in these structures in recent experiments, 
and clarify the degree of validity of the method of moment in an inhomogeneous system. 
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INTRODUCTION 



Nuclear spin polarization and dynamica** in semi- 
conductor nanostructures such as quantum wells and 
quantum dots have attracted increasing attention in re- 
cent years. For example, electrical transport experi- 
ments have demonstrated dynamical nuclear spin polar- 
ization near tunnel junctions, quantum point contacts, 
and coupled quantum dots^s^i^ Optical pumping nu- 
clear magnetic resonance (NMR) technique has been used 
to explore the local electronic state in 2D electron gas 
in the quantum Hall regime by measuring the Knight 
shift and the relaxation time Ti^S, Nuclear spin diffu- 
sion has been found to play an important role in the 
heat capacity anomaly^ at filling factor v = 1, which 
may have originated from a Skyrme solid-liquid phase 
transition. Time-resolved optical measurements in mag- 
netic and non-magnetic semiconductor heterostructures 
also clearly demonstrate strong influences of nuclear spins 
on the confined electron spin dynamics ^Siiii^ 

Both nucleaBi^ii^ii^ and electron spini^ in semiconduc- 
tors have been proposed as the potential quantum bit 
for quantum computing architectures, and nuclear spins 
also are suggested as quantum memoryiii At low tem- 
peratures, the hyperfine interaction between electron and 
nuclear spins could be the dominant decoherence mech- 
anism for both types of spins Because of the con- 
fined nature of electrons in such devices, the hyperfine 
coupling acquires a strongly local characteristics. To 
achieve detailed understanding of electron and nuclear 
spin coherence, a careful study of nuclear spin dynamics 
in these semiconductor heterostructures is imperative. 

One of the nuclear spin relaxation channels is spin dif- 
fusion, which reduces local nuclear polarization through 
direct or mediated spin-spin interaction. Nuclear spin 
diffusion (NSD) was first introduced by Bloembergen 
to explain the measurements of spin-lattice relaxation 
time Ti in ionic crystals in the presence of paramagnetic 
impurities. He suggested that NSD could be induced 
by the mutual nuclear spin flip-flops through magnetic 
dipole-dipole interaction among nuclear spins. Since then 
many calculations^-^^ have been made for the NSD coef- 



ficients. Similar results were obtained via a variety of ap- 
proaches, as these calculations all deal with pure dipole- 
dipole interactions. 



In this paper we present detailed calculations of NSD 
coefficients in semiconductor quantum dots. Although 
the formulation is general, we will concentrate on GaAs 
based dots and wells which are of great experimental in- 
terests. Direct measurements of the NSD coefficients has 
been done using optically pumped NMR for bulk GaAs 
and AlGaAs. It was estimated that the NSD coefficient 
in bulk GaAs is in the order of 10~^^cm^/s for the ar- 
senic nuclei24 and ^ 10~^'^cm^ /s for nuclei in the Al- 
GaAs barrier Our objective in the present study is not 
to accurately predict the numerical values of the NSD co- 
efficients in the nanostructures. Instead, we would like to 
assess how they are modified compared to the bulk ma- 
terials. Specifically, our present focus is on how the hy- 
perfine interactions affect the diffusion coefficients, since 
the confined electrons in these materials have nonuni- 
form wavefunctions, which lead to nonuniform coupling 
to the nuclear spins through the Fermi contact interac- 
tion. Since the hyperfine interaction is much stronger 
than nuclear dipole-dipole interaction wherever the elec- 
tron wavefunction is not negligible, we expect that nu- 
clear spin diffusion could be strongly affected. 



The paper is organized as follows. In section II we 
introduce the moment method^^ and the density matrix 
methodr^iSi which we use to calculate the nuclear spin 
diffusion coefficients. We then discuss how to adapt these 
methods to the inhomogeneous situations of quantum 
dots. In section III, we give numerical results from both 
methods and compare the two approximations. We also 
explore the experimental relevance of our results. Finally 
some further discussion and conclusion are presented in 
section IV. 
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II. FORMULATION 
A. Moment method 

In our calculation, we assume a finite static magnetic 
field Bq along the z direction.^* Under this condition non- 
secular terms of dipolar Hamiltonian can be dropped due 
to consideration of energy conservation, so that the direct 
magnetic dipolar Hamiltonian can be written as^S 

Hi = -7/nBo^/.,+^B,, (2/,,/,, -/,+/,_),(!) 

1 



(2) 



Here 77 is the gyromagnetic ratio of nuclear spin /, Rij is 
the distance between two nuclei located at positions R; 
and Rj , 0ij is the angle between Ry and z direction, and 
'^i=ij stands for the summation over all the spin pairs 
except i — j. We will only consider the dipolar coupling 
among the same nuclear species. Effects of different spin 
species will be briefly discussed in section IV. 

The moment method was designed to study linear re- 
sponse of the spin systemSS, such as the susceptibility of 
the nuclear spin system, by applying a small space- and 
time-dependent magnetic field: 

b{x, t) = BiCoa{Lut)s'm{qx) , 

so that the response of the spin system can be evaluated. 
To study spin diffusion, the perturbing field is along the 
same direction as the static magnetic field. The perturb- 
ing Hamiltonian then takes the form 



-jihBicos{ujt) sin(ga;i 



The 2nth moment is defined asS^t^ 



Mo 



(3) 



(4) 



where a and b are the eigenstates of the unperturbed nu- 
clear spin Hamiltonian and Ea and Eb are the associated 
eigenvalues. The moments contain information on the 
shape of the resonance absorption curve for the whole 
ensemble of nuclear spins. A common practice is to as- 
sume a particular line shape with some unknown param- 
eters, then calculate the first few moments to determine 
these parameters. ■^^ In general, the calculation of M2n 
is rather complicated. However, knowing the first two 
moments is usually enough to determine the line shape 
approximately. In the present situation, after substitut- 
ing Eq. ||2Jl into Eq. ^ , we obtain 



q 



Tr{[i/,/,,][i?, 



2h? 

q 



E.Tr{4} 



(5) 



2^4 



E,Tr{4} 



where Tr represents the thermal average of the operators, 
— x'^ is the difference of the Cartesian coordi- 



nates at nuclear sites H 



and R, , and Greek letters stand 



for the X, y and z directions. In deriving Eq. (|SJ) and 
Eq. ((HJ, it is assumed that the nuclear spins are macro- 
scopically homogeneous so that liz commutes with 
the total Hamiltonian. 

NSD coefficients can be calculated starting from the 
general spin diffusion equation. 



9M(r,t) 
di 



dxt^dx'^ 



(7) 



The diffusion of nuclear magnetization occurs as a result 
of a spatially inhomogeneous initial condition of the mag- 
netization. As we mentioned above the physical mecha- 
nism of NSD is the nuclear spin flip-flops. For a known 
line shape, we can calculate all the moments and evalu- 
ate the spin diffusion coefficients. In most cases the line 
shape can be approximated with a Gaussian. Using the 
Fourier transformed diffusion equation = Dq^, where 
T is the polarization relaxation time; the spin-diffusion 
coefficient D can be expressed in terms of M2 and M 
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If M4^^/3(M 



MM\2 



'Ml; 



(8) 



2 7 is much greater than 1 (corresponding 



to a long tail for the absorption line shape), the Gaus- 
sian approximation becomes inappropriate. A truncated 
Lorentzian shape with a large cutoff frequency is usually 
assumed in such a situation. The spin diffusion coeffi- 
cient D is now 



— 



2\/3 q^ 



(9) 



Since both M2 and M4 are proportional to g^, the diffu- 
sion coefficients in expressions © and ® are indepen- 
dent of q. Notice that the two approximations of line 
shape lead to almost the same numerical results for nu- 
clear spin diffusion coefficients, thus we adopt the Gaus- 
sian line shape (Eq. (H)) throughout this study. 

In the present study we apply the moment method 
to study nuclear spin diffusion in a quantum dot where 
trapped electrons are confined in all three dimensions. 
For simplicity we assume that there is only one electron 
in the dot. The nuclei-electron hyperfine interaction is 
given by 

Bh = ^A(Ri)Ii.S, (10) 

i 

A(RO = ^jaeh'l^mi'- (11) 

Here 7e is the gyromagnetic ratio of the electron in the 
dot, and ^' is the electron wavefunction. In Eq. H1(J|) 
we have ignored the nuclei-electron dipolar interaction. 
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which is much weaker than the contact hyperfine inter- 
action for any finite electron-nucleus distance. 

The nuclear Zeeman energy splitting is about 0.2 
percent of the electron Zeeman energy in GaAs quan- 
tum dots. Furthermore, the electrons in a quantum 
dot has discrete energy spectrum. There is no small 
change of electron kinetic energy that could facilitate 
spin-dependent scattering. Thus direct spin flip-flops be- 
tween the electron and nuclei are largely suppressed in 
strong magnetic fields due to violation of energy conser- 
vation. Here we also neglect any phonon effect since it 
involves a higher order process and is not essential in 
the low temperature limit. The hyperfine interaction in 
Eq. (|10|l can now be reduced to the following effective 
Hamiltonian (assuming electron spin is fully polarized. 
A reduced electron spin polarization will uniformly re- 
duce the strength of Hh): 



(12) 



where and Ai — v4(Ri), and the total Hamiltonian of the 
nuclear spin system is 



Hm = Hi + Hh, 



(13) 



where Hj is the nuclear spin Hamiltonian given in 
Eq. 1^. We notice here that similar approximation of 
neglecting electron-nucleus spin flip-flop has also been 
used to calculate the electron spin spectral diffusion in- 
duced by nuclear spin flip-flops.^" In Eq. p2|l we have 
ignored the spin dynamics of electron, and assumed that 
the electron has been fully polarized. Even if the aver- 
age electron polarization is zero, the calculation of the 
fourth moment in Eq. © would still be non-vanishing, 
since the trace in Eq. (|SJ| involves terms like (S*^) = 1/4 
and (Si) = 1/16. 

The calculation of moments has to be modified in the 
case of a quantum dot. In a homogeneous nuclear spin 
system, the sum over nuclear spin site index i in Eqs. Q 
and jnil is trivial because it means calculating the average 
over the whole homogeneous sample. For the inhomoge- 
neous system considered in the current study, we approx- 
imate the sum over i with the method of coarse graining 
where the sum is evaluated over a few neighboring lat- 
tice points. Such coarse graining is justified since the 
strength of magnetic dipolar interaction decreases quite 
rapidly (l/r^). 

The calculation of the moments is greatly simplified 
at the high temperature limit fc^T 3> h'yjBo, which ap- 
plies in most low temperature experiments (^ 100 mK 
electron temperature), since the nuclear Zeeman energy 
is at the order of 1 mK/Tesla. At the high temperature 
limit we can neglect the Boltzmann factor in the thermal 
averages. The actual evaluation of the commutators and 



traces is long but straightforward. The final results are: 

3 'J 



Tr{[i/M, [HMjrMHm. [HmJj.]]} 



Tr{4} 



= M, 



DD 



A4(15) 



where 
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-t- 4:BijBikBjk{2Bij — Bik — Bjk)} —I (/ + 1) 



/(/ + 1)(16/2 + 16/-7)B4 



IiI + l)Bf^iA-A,r^ 



(16) 
(17) 



Here ^^.j-j ^-j means summation of k over all the lattice 
points except i and j. Mjjjj and Af/j are the dipole-dipole 
contribution and hyperfine contribution to the fourth mo- 
ment, respectively. Our results agree with Redfield and 
Wu's resultsSSi if we set A to be zero. 



B. Density matrix method 

Since the moment method is designed for study of ho- 
mogeneous bulk system, it is important to corroborate 
our results on inhomogeneous systems with a different 
approximation. As a comparison, we calculate diffusion 
coefficients using the density matrix method, ^^'^^ which 
is more straightforward in terms of its physical picture. 
We assume that the density matrix of nuclear spin sys- 
tem can be expanded in terms of a group of orthogonal 
operators {i — 1, N), where N is the number of nuclear 
spins in the system: 



P 



and 



Tr{/,,/,J = %Tr{4}. 



(18) 



(19) 



This choice for the nuclear spin density matrix is a good 
approximation at the high temperature limit, which is 
usually satisfied by the systems we are interested in. The 
total Hamiltonian for the nuclear spin system is 

-ffoM = Ho + Hi, 

Ho = - lihBo)U, + 2 B,,I,Jj„ 



Hi - -^Bi-jli+Ij 



(20) 



Here we take the nuclear flip-flop term (which accounts 
for spin diffusion) as a perturbation. 
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The nuclear spin density matrix can be conveniently 
calculated in the interaction picture 



Pit) = e"°'p{t) 



-iHot 



(21) 



We have used " ^ " to represent the operators in the in- 
teraction picture. The temporal dynamics of the density 
matrix in the interaction picture is governed by the flip- 
flop term in the full Hamiltonian 



A second-order calculation leads to 



(22) 



~p{t)=i[p{t),Hi{t)]+i^ / dT[H^{t),[H^{t-T), pirn- 
Jo 

(23) 

Substituting Eq. UHl) into Eq. |(23l, we find 

afc(t) (24) 



with 



-iTr 



[Hiit),h,]Ikz 



drTr 



[Hi{t),hz][Hi{t-T),U, 



.(25) 



One can easily show that Wki = Wik by noting that 
the trace is invariant under the cyclic reordering of op- 
erators. Wik describes the flip-flop rate of two nuclear 
spins at site i and k. Substituting the density matrix 



p{t) ~ o.i{t)hz into the equation of motion of the 
local nuclear magnetization 



d_ 

dt 



ihz) = Tr{p(t)4,} = Tr{^(0/fe^} 



(26) 



and performing Taylor expansion around the space point 
of the fcth nucleus, we find 



(27) 



In writing Eq. I|26|l we have used the commutation re- 
lation [Ikz,Ho] — 0. It is easy to show that Wkk ~ 0, 
because it involves a summation of a fast oscillatory func- 
tion that averages to zero over many nuclear sites. Phys- 
ically, Wkk corresponds to energetically impossible pro- 
cesses and has no physical meaning. It then follows that 
the zero-order term in the Taylor expansion does not con- 
tribute to spin diffusion. The first-order term also van- 
ishes because of the crystal symmetry,?" 

To calculate the diffusion coefficients we need to find 
the flip-flop rates Wik- The explicit calculations of the 
traces for an arbitrary nuclear spin in Eq. H25|) are quite 
complicated. In the following we consider the particular 



situation of spin 3/2 nuclei, which is the case for GaAs 
quantum dot. Calculating the trace for I — 3/2, we 
obtain 

TT[[Hi{t),Ikz][Hi{t-T),I,,]} = 2i?2,cos(^T) 

[cos{2Bik,nT) + COs{6BikmT)] , (28) 

m(i,k) 

where f{x) = 34-1- 48cos(a;) -I- 18cos(2x). Here we have 
used the deflnition Aik = A^-Ak and B.^km = Bi„^-Bum- 
Finally we get the expression of Wik 



B 



2 



,A 



W,k = ^ dT cosier) f{AB,kT) 

X Jl C0s{4:BikmT) COs{2BikmT) , (29) 

m{i,k) 



which would then allow us to calculate the NSD coeffi- 
cient of the system. 



III. NUMERICAL RESULTS 

A. NSD in bulk system 

Before presenting our numerical results for a quantum 
dot, we first estimate the NSD coefficients for pure nu- 
clear spin dipole-dipole interaction using Eq. with 
the moment method. Notice that the hyperfine inter- 
action does not change the first moment. The sum- 
mations in Eq. (|16|) can be easily done over the nuclei 
in a face-centered cubic structure (for GaAs). Since 
the dipole interaction decays as r""^, the summations 
converge quite rapidly. A numerical calculation yields 
D-\= 0.29jjh/aG.As and D^^ = 0.167|ft/aGaAs for 
/ = |, where the lattice constant acaAs — 5.65 A. These 
values are comparable to Lowe and Cade's result^ for 
spin one half in a simple cubic structure. For the spe- 
cific example of ^^As nuclei, where 7/ = 4.58 x 10^^, 
= 1.1 X lO-i^cmVs and D^"^ = 6.3 x IQ-^cm^/s. 
Spin diffusion is faster along the z direction because the 
dipolar interaction is stronger along the external mag- 
netic field direction according to Eq. (0) . Specifically, the 
dipolar coupling coefficient is proportional to the magni- 
tude of 1 — 3cos^6'jj. Along z direction this value is -2, 
while it is 1 along x or y direction. In the following dis- 
cussion we use Dq^ to represent the NSD coefficient for 
pure dipole-dipole interaction in the absence of inhomo- 
geneity. 

To calculate the NSD coefficients with the density ma- 
trix method, we have to evaluate the integral in Eq. H29|l . 
This can be done by first changing the upper limit of the 
integration to infinity because the integrand is a product 
of many cosine functions that has a sharp spectral peak 
near t = 0, so that changing the integration upper limit 
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only introduces a negligible error ifi* We thus have 



N 



du ]^cos(aiu) 

1=1 

j exp Mn cos(aiM) 



du exp — au \ = —\ — , 



(30) 



where a = af. In the second step of the calcula- 
tion in Eq. H3l)|l we have expanded the integrand around 
u = and kept only the terms to the order 0{v?). This 
approximation is in the same spirit as the steepest de- 
scent method. Using this approximation we find Eq. 1291 
takes the following form 



(0) 



,-7,. 



(1) 



(2) 



F, 



ik 



17V27r 
5 



Biki^fk + dik)' 



(2) 



5 

/2^ 



i^ki^. 



Qik 



80 



10 

E 



^ik i^ik 

(Bip — 



256BI 



9ik} 



9ik) 



Bk, 



kp) 



(31) 



The calculated Wik can then be inserted into Eq. H27|) 
to obtain the diffusion coefficients. For pure dipolar 
interaction we find D^^ — 0.497|/i/aGaAs and I?^^ — 
0.2l7jfi/aGaAs- These calculated D^^ and D^'^ are nearly 
twice as large as the results given by the moment calcula- 
tions, although we do find that D^^ is greater than D^^, 
similar to the results of the moment calculations. 



B. NSD in a quantum dot 

We now include electron-nuclear spin hyperfine inter- 
action in our calculation of the NSD coefficients. To 
study the effects of hyperfine interaction, we need knowl- 
edge of the electronic wavef unctions. The ground state 
electron wavefunction in a 2D gated GaAs quantum dot 
can be approximated by 




(32) 



where zq is the quantum dot thickness, Iq is the 
Fock-Darwin radius, and rp is the electrostatic lateral 
parabolic confinement radius. The value of the F-point 
Bloch function tt(r) at nuclear sites can be deduced from 
experimental measurementst^i 

The calculated NSD coefficients using the moment 
method and density matrix approach share several com- 
mon characteristics. Figure ^ shows both D^^ / Dq^ and 




z(10 m) 



z(10 m) 



FIG. 1: The ratio of NSD coefficients D"" /D^" as a function 
of spatial coordinate z for various Fock-Darwin radii, where z 
is the perpendicular distance from the center of the dot. The 
left panel (a) shows the results using moment method while 
the right panel (b) represents those obtained with density ma- 
trix method. In all these calculations we assume a quantum 
dot with thickness zo = 10 nm. The solid line (upward trian- 
gle) represents Zo = 30 nm for 7?^^ (D^^) and the dashed line 
(downward triangle) describes Iq — 80 nm for (D^^). 



jjxx I jjxx ^ functions of spatial coordinate (along the 
external magnetic field) z for two different quantum dot 
thickness. In Fig. |3 we plot D'^^/D^^ and 
as functions of the radial displacement r (perpendicu- 
lar to the external magnetic direction) for different Fock- 
Darwin radius Iq. The curves in both figures show similar 
behaviors. The suppression of spin diffusion due to hy- 
perfine interaction at or near the center of the quantum 
dot could be so significant that D^^ and D^^ is only a 
few percent of Dg^ and Dg^. Figures ^ and |21 also show 
that with both methods the suppression of spin diffusion 
decreases as the dot size becomes larger, which can be 
explained by noting that hyperfine interaction strength 
decreases for larger dots. To further illustrate this point, 
in Fig. O we show the diffusion coefhcient D^^ at the 
center of the quantum dot as a function of Fock-Darwin 
radius /q- Similar results (which are not shown in Fig.|31) 
are found for D^^ as well. An additional feature of Fig. 
is that, at the boundary of the dot along z direction, the 
NSD coefficients increase to Dq rapidly. This behavior 
is due to our assumption that the electron wavefunction 
outside the dot is zero. It is also noticed that 13^^ has 
a stronger suppression than D^^ , which is illustrated in 
both Fig. n and Fig. [fusing both calculation methods. 

There are several interesting differing features in the 
results of the two methods in addition to the different 
magnitudes of the diffusion coefhcients given by the two 
methods as we have discussed previously for homoge- 
neous systems. Firstly, moment method typically leads 
to stronger suppressions at the center of the dot. Sec- 
ondly, Fig. n particularly panel (b), indicates that the 
suppression of NSD is not the strongest at the center 
of the quantum dot. Instead it decreases from the cen- 
ter (albeit only slightly in some cases), reaches its mini- 
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FIG. 2: The ratio of NSD coefficients D"" /i?^'" as a function 
of spatial coordinate r for various quantum dot thickness, 
where r is the radial displacement in the 2D plane. The origin 
of the coordinate system is located at the center of the dots. 
The left panel (a) are the results of moment method and the 
right panel (b) shows those from density matrix method. In 
all these calculations we have used a quantum dot with Fock- 
Darwin radius Iq = 50 nm. The solid line (upward triangle) 
describes zo = 5 nm for D^^ (D^^) while the dashed line 
(downward triangle) shows zq — 15 nm for D^^ (D^^). 




FIG. 4: The ratio of NSD coefficients D^^/Dg^ (calculated 
using the density matrix method) as a function of spatial co- 
ordinate r for strongly confined dots (panel a) and quantum 
wells (panel b) with Gaussian shaped carrier spin polarization 
(see discussion in the text). In the case of quantum dots we 
assume zq = 10 nm and Iq — 25 nm. In the calculation of 
quantum well with thickness of 15 nm, we assume the effec- 
tive radius of the Gaussian shaped carrier spin polarization is 
12 ^m. In both panels the thinner lines show the variations 
of hyperfine coupling strength as a function of r. In panel b, 
the two solid and dashed lines represent different polarization 
magnitude. 



mum at an intermediate position, then starts to rise again 
near the edge of the dot. This feature is stronger in 
the results obtained through the density matrix method 
than those from the moment method, and is stronger for 
the z-direction diffusion (Fig. ^) than in-plane directions 
(Fig. [SJ. In fact the characteristic exists, although only 
weakly, in Fig. [3b) for D^^ while it is not present in 
Fig. Eta). 

The presence of off-center local minimum in nuclear 
spin diffusion in a quantum dot is physical. The sup- 
pression of spin diffusion is determined by the difference 
of the inhomogeneous hyperfine coupling at two nuclear 
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FIG. 3: The ratio of NSD coefficients D'"'' / Dq'' at the center 
of the dot as a function of Fock-Darwin radius lo with three 
difi'erent dot thickness. The results are obtained with the 
density matrix method. 



sites cx - IvE'jp V|*ip • (r^ - r^) instead of the 
coupling constant alone (see Eq. (|17|l and Eq. (|29|l ). In 
essence, in the flip-flop processes that account for nuclear 
spin diffusion, energy must be conserved. If the hyper- 
fine coupling strengths are not the same at the two lat- 
tice sites, the effective Zeeman energies of the two nuclear 
spins are different, so that extra energy must be absorbed 
or emitted (from the overall dipolar energy reservoir, for 
example) to compensate for the difference. Apparently 
small energy differences should result in greater flip-flop 
rates, hence faster diffusion. This is the basic reason why 
spin diffusion is suppressed in an inhomogeneous system. 
Since the energy difference is proportional to both the 
magnitude and the gradient of the electron wavefunction, 
the strongest suppression of NSD coefficient could occur 
either at the center of a quantum dot or its "waist" , where 
the gradient is the largest (notice that the nuclear spins 
are on discrete sites, thus near the center of the quantum 
dot the hyperfine energy difference is generally nonvan- 
ishing, for a very small dot it could even be a maximum 
depending on the form of the electron wavefunction) . For 
an electron confined in a quantum dot geometry, the en- 
velope wavefunction is highly non-uniform and can be 
approximated with Ea.f Along the z direction the 
changes of this hyperfine coupling between neighboring 
nuclear sites are quite large. This is the reason that there 
is a sharp minimum of the diffusion coefficient as a func- 
tion of z at a nonzero z. On the other hand the confine- 
ment is not that strong in the r direction especially for 
larger dots, so the feature is not as obvious. 

To illustrate the previous discussion in more detail we 
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consider a small quantum dot with = 25 nm and zq = 
10 nm. Figure EJa) clearly shows a sharp olf-center local 
minimum of the NSD coefficient D^^, corresponding to a 
strong suppression of diffusion, which is weak in Fig. |21 
Quite interestingly a ring structure of nuclear spin polar- 
ization has been observed in ferromagnet /semiconductor 
heterostructures by spatially modulating the excitation 
intensity. '^^ The experiment uses a laser pulse with a 
Gaussian cross section, which we believe is key to the 
ring structure. The inhomogeneous power input induces 
non-uniform carrier polarization strength at the inter- 
face which in turn leads to inhomogeneous hyperfine cou- 
plings. As we have discussed, the suppression of diffusion 
would be the strongest at some position between the cen- 
ter and the boundary. In Fig.^fb), we show the suppres- 
sions of NSD coefficient D^^ for two different polarization 
magnitudes in the micrometer size. The results are quite 
similar to the strongly confined quantum dots. At low 
temperatures an important spin relaxation mechanism 
is spin diffusion. The existence of the local minimum 
of diffusion coefficients contributes to a maximum of nu- 
clear spin polarization since the diffusion is slowest at the 
point. In other words a ring structure could very well be 
present, as what was observed experimentally. 

IV. DISCUSSION AND CONCLUSION 

In the present study we have investigated the dipole- 
dipole interaction among like nuclear species (Eq. l|T]l). 
Interaction between unlike nuclear spins have been ne- 
glected. Under the assumption that the magnetic field 
is not weak, this should be a good approximation. How- 
ever there is the so-called indirect interaction (RKKY)^ 
in highly disordered samples where spin-flip scattering 
has measurable physical effects. In this regard, the cou- 
pling between different nuclear species may have non- 
negligible effects. It should be mentioned that M2 and 
M4 in Eq. H14|l and Eq. (|15|l do not change without direct 
spin interaction. However, the evaluation of M4 becomes 
extremely complicated if the indirect coupling between 
unlike spin species is included. We did not study this 
aspect of spin diffusion in the current paper. 

We have calculated NSD coefficients for Arsenic nuclei 
in this paper. GaAs has a zinc blende structure with 
50% ^^As (which is the only stable Arsenic isotope). In 
natural GaAs samples, there are two isotopes of Gallium, 
^iGa (19.8%, 7/ = 8.16 x lO^l/s • G) and ^SGa (30.2%, 
7/ = 6.42 X lO'^l/s • G). In the barrier region, the Ga 
concentration is even lower with the introduction of 10% 
to 15% of Al in place of Ga. An evaluation of the NSD 
coefficients for Ga would have to account for the random 
distribution of different Ga isotopes on the fee lattice. 
Here our emphasis is the effect of inhomogeneous hyper- 
fine interaction on NSD. Furthermore, this nonuniform 
hyperfine coupling, in the form of Eq. 1)12(1 , cannot com- 
pensate for the difference in Zeeman energy of different 



nuclear species, so that the inter-species NSD is unlikely. 
For example, the effective hyperfine magnetic field seen 
by nuclei at the center of a quantum dot is only a few 
tens Gauss, which is usually much less than the external 
field. Thus the inter-species NSD is basically impossible 
in a finite magnetic field, and we do not have to consider 
the Ga nuclei when studying NSD of the As nuclei. 

Recently spin diffusion suppression by nonuniform field 
has been found^ in silica samples where an inhomoge- 
neous magnetic field was generated by a ferromagnetic tip 
of a magnetic resonance force microscope. It was found 
that spin relaxation rates T^^ are significantly reduced 
due to the suppression of nuclear spin flip-flop processes. 
In solids with paramagnetic impurities, inhomogeneous 
internal field could also be generated by dipole-dipole in- 
teraction between the impurity and its neighboring nu- 
clear spins, in which case a barrier to NSD can also be 
formed. 

In this study we have presented two methods to study 
the suppression of spin diffusion. In the moment method 
originally developed for homogeneous bulk material;^ a 
small spatially and temporally varying perturbation is 
added to the total Hamiltonian to generate linear re- 
sponse. The spatial variation of the perturbation is as- 
sumed to be smooth compared to electron wavefunction 
variation (long wavelength approximation). To apply it 
to nanostructures like quantum dots, we combine it with 
a coarse-graining approximation. In this method we have 
to assume the sizes of quantum dots being considered 
are relatively large. In principle the moment method is 
not designed for the calculation of strong spatial varia- 
tion of diffusion. This partly explains the less prominent 
local minimum of diffusion coefficients as a function of 
spatial coordinates. On the other hand the density ma- 
trix method is more straightforward and keeps more local 
features in the evaluations. With our density matrix cal- 
culation no assumption of any particular line-shape is 
necessary. We simply start from the full Hamiltonian, 
and then use the second-order perturbation to find the 
evolution of the local nuclear magnetization. Neverthe- 
less the strong suppressions of NSD coefficients for small 
dots appear in both calculations. 

To conclude we have presented detailed study of nu- 
clear spin diffusion under the influence of inhomogeneous 
contact hyperfine interactions in GaAs based nanostruc- 
tures. Our results show that there are strong suppres- 
sions of spin diffusion at the center and the waist of 
a quantum dot or quantum well where electron prob- 
ability and/or gradient of electron probability is large, 
which is consistent with experimental observations in 
such structures. ^'^^ The numerical results given in Sec- 
tion HI show that NSD coefficients could be suppressed to 
as small as a few percent of Dq^ . Our results clearly show 
that non-uniform electron distribution can help maintain 
desired nuclear spin polarization in these semiconductor 
nanostructures. 
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